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1.1 Methods of Describing Fluid Motion 

To describe the fluid motion, one needs to know the variations of physical quantities such 
as , density, velocity, pressure, temperature, stresses, etc., as functions of time, everywhere 
within a certain spatial region. There are two ways to describe the fluid motion. One is 
called Lagrangian, where one follows all fluid particles and describes the variations around 
each fluid particle along its trajectory. The other is Eulerian, where the variations are 
described at all fixed stations as a function of time. In the second, different particles pass 
the same station at different times. 

Consider first the Lagrangian description. 

1.1.1 Lagrangian description 

A particle is identified by its initial position at time to, 

x = (x ,y ,z ) (1-1-1) 

Let 

x = x(x ,t) (1.1.2) 

be the position of the same particle at the later time t. By definition, at t — to, x (xq, to) — %o- 
Instead of the initial position xq, one can use any three quantities a, b, c, which are uniquely 
related to x : 

x = x (a) a = (a, b, c) (1.1.3) 
to identify the particle. Thus the path of a particle identified by a is given by : 

x = x(a,t). (1.1.4) 

Note that a, t are the independent variables, and x = (x, y, z) are dependent variables. 
From the particle position we can calculate the particle velocity: 



dx 



(1.1.5) 
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as well as the particle acceleration: 



dt 



d 2 x 



dt< 



(1.1.6) 



Other physical quantities such as density p and pressure p can also be expressed in terms of 
a, t, e.g., 

p = p(a,t), p = p(a,t). (1-1-7) 
1.1.2 Eulerian Description 

Here we specify the fluid properties (density, velocity, pressure...) at a fixed point and a 
chosen time. Thus we are dealing with the time records of a property at a fixed measuring 
probe. 

q(x,t) , p(x,t) , p(x,t) , etc. (1.1.8) 

Once a physical property of all particles are known at all times, how do we get the property 
at all fixed points at all times? In other words, how are the Eulerian and Lagrangian systems 
related? Suppose first that the Lagrangian information of the velocity of all particles are 
known : 

q = q(a, t) 



dx 



= q(a,t) 



dt 

x — x (a, t) 



then we can integrate the definition 

for to get the position of a particle, 
so that 

x(a, 0) = a. 

The result can be inverted in principle to get 

a — a(x, t) 
as long as the Jacobian of transformation 



1.1.9) 



(1.1.10) 



'i.i. ir 



(1.1.12) 
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7^0. 



1.1.13) 



does not vanish. Once the position of a particle at time t is known, any physical property 
of a particle can be translated to the property at certain position in space. For example, by 
substituting (1.1.12) into (1.1.9), we get 



q = q(x, t) , 



1.1.14) 
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which gives the Eulerian velocity, i.e., the velocity field at all fixed points. Similarly from 
the Lagrangian pressure p((a,t) we get the Eulerian pressure 

p((a(x,t),t) =p(x,t) 

Consider next that the Eulerian velocity field is known everywhere: 

q = q(x,t) (1.1.15) 

Since at t the particle a is at x (a, t); the Lagrangian velocity of particle a is the same as the 
Eulerian velocity. Therefore, 

Ox 

q(x,t) = —, (1.1.16) 

This equation is a system of nonlinear ordinary differential equations. If they can be solved 
for x = x (a, t) subjected to the initial condition x = xo (a) , t = t , we then have 

Ox 

— = q[x(a,t),t] 

= <lL( a ,t) = ^Lagrangian 

For the rate of variations we need to calculate derivatives of some physical quantity 
A (x, t) whose Eulerian information is given. The rate of variation following a fluid particle 

dA = <M , dA 8A , dA 

dt s dt ^ dx t dt s dt ^ dxi (1.1.17) 

(Lagrangian) (Eulerian) 



The operator 

D , 

Di = m +q 'd7, (LL18) 

has different names: the substantial, or total, or material derivative. The second term is 
referred to as the convective rate of variation since it arises because the fluid particle moves 
to new places. In particular, the acceleration of a fluid particle is 

A more physical interpretation of Df /Dt can be derived by using the Eulerian picture 
directly. Consider f(x,t). After 5t, f(x,t) is changed to f(x + q5t,t + St). Therefore, the 
total change of / is 

5f = f[x + q5t,t + 5t]-f(x,t) 

= f(x,t)+(q-Vf + ^)st-f(x,t) = ^5t. 



It follows that 

s_i 

5t 



= w = l +? ' v/ (112 °) 



where |£ is the instantaneous change, and q ■ V/ is the convective change. 



